An integrated metagenomic, metabolomic and transcriptomic survey of Populus across genotypes and environments.

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

Merit of integrating in situ transcriptomics and anatomical information for cell annotation and lineage construction in single-cell analyses of Populus.

Cell type annotation and lineage construction are two of the most critical tasks conducted in the analyses of single-cell RNA sequencing (scRNA-seq). Four recent scRNA-seq studies of differentiating xylem propose four models on differentiating xylem development in Populus. The differences are mostly caused by the use of different strategies for cell type annotation and subsequent lineage interpretation. Here, we emphasize the necessity of using in situ transcriptomes and anatomical information to construct the most plausible xylem development model.

Genome-Wide Identification and Analysis of the EIN3/EIL Transcription Factor Gene Family in Doubled Haploid (DH) Poplar.

Ethylene (ET) is an important phytohormone that regulates plant growth, development and stress responses. The ethylene-insensitive3/ethylene-insensitive3-like (EIN3/EIL) transcription factor family, as a key regulator of the ET signal transduction pathway, plays an important role in regulating the expression of ET-responsive genes. Although studies of EIN3/EIL family members have been completed in many species, their role in doubled haploid (DH) poplar derived from another culture of diploid Populus simonii × P. nigra (donor tree, DT) remains ambiguous. In this study, a total of seven EIN3/EIL gene family members in the DH poplar genome were identified. Basic physical and chemical property analyses of these genes were performed, and these proteins were predicted to be localized to the nucleus. According to the phylogenetic relationship, EIN3/EIL genes were divided into two groups, and the genes in the same group had a similar gene structure and conserved motifs. The expression patterns of EIN3/EIL genes in the apical buds of different DH poplar plants were analyzed based on transcriptome data. At the same time, the expression patterns of PsnEIL1, PsnEIN3, PsnEIL4 and PsnEIL5 genes in different tissues of different DH plants were detected via RT-qPCR, including the apical buds, young leaves, functional leaves, xylem, cambium and roots. The findings presented above indicate notable variations in the expression levels of PsnEIL genes across various tissues of distinct DH plants. Finally, the PsnEIL1 gene was overexpressed in DT, and the transgenic plants showed a dwarf phenotype, indicating that the PsnEIL1 gene was involved in regulating the growth and development of poplar. In this study, the EIN3/EIL gene family of DH poplar was analyzed and functionally characterized, which provides a theoretical basis for the future exploration of the EIN3/EIL gene function.

Impacts of AlaAT3 transgenic poplar on rhizosphere soil chemical properties, enzyme activity, bacterial community, and metabolites under two nitrogen conditions.

Poplar stands as one of the primary afforestation trees globally. We successfully generated transgenic poplar trees characterized by enhanced biomass under identical nutrient conditions, through the overexpression of the pivotal nitrogen assimilation gene, pxAlaAT3. An environmental risk assessment was conducted for investigate the potential changes in rhizosphere soil associated with these overexpressing lines (OL). The results show that acid phosphatase activity was significantly altered under ammonium in OL compared to the wild-type control (WT), and a similar difference was observed for protease under nitrate. 16SrDNA sequencing indicated no significant divergence in rhizosphere soil microbial community diversity between WT and OL. Metabolomics analysis revealed that the OL caused minimal alterations in the metabolites of the rhizosphere soil, posing no potential harm to the environment. With these findings in mind, we anticipate that overexpressed plants will not adversely impact the surrounding soil environment.

Genome-Wide Identification and Expression Profiling of Velvet Complex Transcription Factors in Populus alba × Populus glandulosa.

The discovery of new genes with novel functions is a major driver of adaptive evolutionary innovation in plants. Especially in woody plants, due to genome expansion, new genes evolve to regulate the processes of growth and development. In this study, we characterized the unique VeA transcription factor family in Populus alba × Populus glandulosa, which is associated with secondary metabolism. Twenty VeA genes were characterized systematically on their phylogeny, genomic distribution, gene structure and conserved motif, promoter binding site, and expression profiling. Furthermore, through ChIP-qPCR, Y1H, and effector-reporter assays, it was demonstrated that PagMYB128 directly regulated PagVeA3 to influence the biosynthesis of secondary metabolites. These results provide a basis for further elucidating the function of VeAs gene in poplar and its genetic regulation mechanism.

Nanoplastic toxicity induces metabolic shifts in Populus × euramericana cv. '74/76' revealed by multi-omics analysis.

There is increasing global concern regarding the pervasive issue of plastic pollution. We investigated the response of Populus × euramericana cv. '74/76' to nanoplastic toxicity via phenotypic, microanatomical, physiological, transcriptomic, and metabolomic approaches. Polystyrene nanoplastics (PS-NPs) were distributed throughout the test plants after the application of PS-NPs. Nanoplastics principally accumulated in the roots; minimal fractions were translocated to the leaves. In leaves, however, PS-NPs easily penetrated membranes and became concentrated in chloroplasts, causing thylakoid disintegration and chlorophyll degradation. Finally, oxidant damage from the influx of PS-NPs led to diminished photosynthesis, stunted growth, and etiolation and/or wilting. By integrating dual-omics data, we found that plants could counteract mild PS-NP-induced oxidative stress through the antioxidant enzyme system without initiating secondary metabolic defense mechanisms. In contrast, severe PS-NP treatments promoted a shift in metabolic pattern from primary metabolism to secondary metabolic defense mechanisms, an effect that was particularly pronounced during the upregulation of flavonoid biosynthesis. Our findings provide a useful framework from which to further clarify the roles of key biochemical pathways in plant responses to nanoplastic toxicity. Our work also supports the development of effective strategies to mitigate the environmental risks of nanoplastics by biologically immobilizing them in contaminated lands.

Integrative omics studies revealed synergistic link between sucrose metabolic isogenes and carbohydrates in poplar roots infected by Fusarium wilt.

Advances in carbohydrate metabolism prompted its essential role in defense priming and sweet immunity during plant-pathogen interactions. Nevertheless, upstream responding enzymes in the sucrose metabolic pathway and associated carbohydrate derivatives underlying fungal pathogen challenges remain to be deciphered in Populus, a model tree species. In silico deduction of genomic features, including phylogenies, exon/intron distributions, cis-regulatory elements, and chromosomal localization, identified 59 enzyme genes (11 families) in the Populus genome. Spatiotemporal expression of the transcriptome and the quantitative real-time PCR revealed a minuscule number of isogenes that were predominantly expressed in roots. Upon the pathogenic Fusarium solani (Fs) exposure, dynamic changes in the transcriptomics atlas and experimental evaluation verified Susy (PtSusy2 and 3), CWI (PtCWI3), VI (PtVI2), HK (PtHK6), FK (PtFK6), and UGPase (PtUGP2) families, displaying promotions in their expressions at 48 and 72 h of post-inoculation (hpi). Using the gas chromatography-mass spectrometry (GC-MS)-based non-targeted metabolomics combined with a high-performance ion chromatography system (HPICS), approximately 307 metabolites (13 categories) were annotated that led to the quantification of 46 carbohydrates, showing marked changes between three compared groups. By contrast, some sugars (e.g., sorbitol, L-arabitol, trehalose, and galacturonic acid) exhibited a higher accumulation at 72 hpi than 0 hpi, while levels of α-lactose and glucose decreased, facilitating them as potential signaling molecules. The systematic overview of multi-omics approaches to dissect the effects of Fs infection provides theoretical cues for understanding defense immunity depending on fine-tuned Suc metabolic gene clusters and synergistically linked carbohydrate pools in trees.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Differences in mineral and osmotic balances enhance zinc translocation in an aquaporin overexpressing poplar.

Zinc (Zn) is an essential micronutrient for plants, but it is toxic beyond a certain threshold. Populus alba (L.) 'Villafranca' clone is known for its good tolerance to high Zn concentration compared to other poplar species. A line of this species overexpressing the tonoplast intrinsic aquaporin AQUA1 gene has showed an improved tolerance to Zn excess in comparison to the wild-type (wt) line. The aims of this work were to: 1) verify if AQUA1 plants can uptake Zn more efficiently after a longer period of exposure; 2) evaluate if a higher Zn uptake in transgenic lines can have negative effects; 3) assess Zn competing elements (iron and manganese), soluble sugars, osmolytes, and potassium to investigate differences in water and osmotic homeostasis between lines. Under Zn excess, AQUA1 plants showed a twofold Zn translocation factor and a higher xylem sap Zn concentration than the wt plants. Transgenic plants preferentially allocated Zn in aerial biomass and this different behaviour matched with modified manganese and iron balances suggesting that the increased Zn uptake might be related to a decrease in iron transport in the transgenic line. Moreover, a higher instantaneous water use efficiency in control conditions and an increase in bark soluble sugars under Zn excess could allow a higher resistance of AQUA1 plants to the water and osmotic perturbations caused by Zn. Indeed, the Zn excess increased the xylem osmolyte content only in wt plants. Further investigations are required to understand the role of AQUA1 in osmotic regulation.

Genome-Wide Characterization and Expression Analyses of Major Latex Protein Gene Family in Populus simonii × P. nigra.

Major latex proteins, or MLPs, are crucial to plants' capacity to grow, develop, and endure biotic and abiotic stresses. The MLP gene family has been found in numerous plants, but little is known about its role in Populus simonii × P. nigra. This study discovered and assessed 43 PtMLP genes that were unevenly dispersed throughout 12 chromosomes in terms of their physicochemical characteristics, gene structure, conserved motifs, and protein localization. Based on their phylogeny and protein structural characteristics, three separate subclasses of PtMLP family were identified. Segmental and tandem duplication were found to be essential variables in the expansion of the PtMLP genes. The involvement of the PtMLP genes in growth and development, as well as in the responses to different hormones and stresses, was demonstrated by cis-regulatory element prediction. The PtMLP genes showed varying expression patterns in various tissues and under different conditions (cold, salt, and drought stress), as demonstrated in RNA-Seq databases, suggesting that PsnMLP may have different functions. Following the further investigation of the genes demonstrating notable variations in expression before and after the application of three stresses, PsnMLP5 was identified as a candidate gene. Subsequent studies revealed that PsnMLP5 could be induced by ABA treatment. This study paves the way for further investigations into the MLP genes' functional mechanisms in response to abiotic stressors, as well as the ways in which they can be utilized in poplar breeding for improved stress tolerance.

Genome-Wide Identification, Expression Analysis, and Subcellular Localization of DET2 Gene Family in Populus yunnanensis.

(1) Background: Brassinosteroids (BRs) are important hormones involved in almost all stages of plant growth and development, and sterol dehydrogenase is a key enzyme involved in BRs biosynthesis. However, the sterol dehydrogenase gene family of Populus yunnanensis Dode (P. yunnanensis) has not been studied. (2) Methods: The PyDET2 (DEETIOLATED2) gene family was identified and analyzed. Three genes were screened based on RNA-seq of the stem tips, and the PyDET2e was further investigated via qRT-PCR (quantitative real-time polymerase chain reaction) and subcellular localization. (3) Results: The 14 DET2 family genes in P. yunnanensis were categorized into four groups, and 10 conserved protein motifs were identified. The gene structure, chromosome distribution, collinearity, and codon preference of all PyDET2 genes in the P. yunnanensis genome were analyzed. The codon preference of this family is towards the A/U ending, which is strongly influenced by natural selection. The PyDET2e gene was expressed at a higher level in September than in July, and it was significantly expressed in stems, stem tips, and leaves. The PyDET2e protein was localized in chloroplasts. (4) Conclusions: The PyDET2e plays an important role in the rapid growth period of P. yunnanensis. This systematic analysis provides a basis for the genome-wide identification of genes related to the brassinolide biosynthesis process in P. yunnanensis, and lays a foundation for the study of the rapid growth mechanism of P. yunnanensis.

N6-methyladenosine mRNA methylation positively regulated the response of poplar to salt stress.

As the most abundant form of methylation modification in messenger RNA (mRNA), the distribution of N6-methyladenosine (m6A) has been preliminarily revealed in herbaceous plants under salt stress, but its function and mechanism in woody plants were still unknown. Here, we showed that global m6A levels increased during poplar response to salt stress. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) revealed that m6A significantly enriched in the coding sequence region and 3'-untranslated regions in poplar, by recognising the conserved motifs, AGACU, GGACA and UGUAG. A large number of differential m6A transcripts have been identified, and some have been proved involving in salt response and plant growth and development. Further combined analysis of MeRIP-seq and RNA-seq revealed that the m6A hypermethylated and enrich in the CDS region preferred to positively regulate expression abundance. Writer inhibitor, 3-deazaneplanocin A treatment increased the sensitivity of poplar to salt stress by reducing mRNA stability to regulate the expression of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Furthermore, we verified that the methyltransferase PagFIP37 plays a positively role in the response of poplar to salt stress, overexpressed lines have stronger salt tolerance, while RNAi lines were more sensitive to salt, which relied on regulating mRNA stability in an m6A manner of salt-responsive transcripts PagMYB48, PagGT2, PagNAC2, PagGPX8 and PagARF2. Collectively, these results revealed the regulatory role of m6A methylation in poplar response to salt stress, and revealed the importance and mechanism of m6A methylation in the response of woody plants to salt stress for the first time.

GWAS supported by computer vision identifies large numbers of candidate regulators of in planta regeneration in Populus trichocarpa.

Plant regeneration is an important dimension of plant propagation and a key step in the production of transgenic plants. However, regeneration capacity varies widely among genotypes and species, the molecular basis of which is largely unknown. Association mapping methods such as genome-wide association studies (GWAS) have long demonstrated abilities to help uncover the genetic basis of trait variation in plants; however, the performance of these methods depends on the accuracy and scale of phenotyping. To enable a large-scale GWAS of in planta callus and shoot regeneration in the model tree Populus, we developed a phenomics workflow involving semantic segmentation to quantify regenerating plant tissues over time. We found that the resulting statistics were of highly non-normal distributions, and thus employed transformations or permutations to avoid violating assumptions of linear models used in GWAS. We report over 200 statistically supported quantitative trait loci (QTLs), with genes encompassing or near to top QTLs including regulators of cell adhesion, stress signaling, and hormone signaling pathways, as well as other diverse functions. Our results encourage models of hormonal signaling during plant regeneration to consider keystone roles of stress-related signaling (e.g. involving jasmonates and salicylic acid), in addition to the auxin and cytokinin pathways commonly considered. The putative regulatory genes and biological processes we identified provide new insights into the biological complexity of plant regeneration, and may serve as new reagents for improving regeneration and transformation of recalcitrant genotypes and species.

Comparative analysis of PLATZ transcription factors in six poplar species and analysis of the role of PtrPLATZ14 in leaf development.

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are a class of plant-specific transcription factor that play a crucial role in plant growth, development, and stress response. However, the evolutionary relationship of the PLATZ gene family across the Populus genus and the biological functions of the PLATZ protein require further investigation. In this study, we identified 133 PLATZ genes from six Populus species belonging to four Populus sections. Synteny analysis of the PLATZ gene family indicated that whole genome duplication events contributed to the expansion of the PLATZ family. Among the nine paralogous pairs, the protein structure of PtrPLATZ14/18 pair exhibited significant differences with others. Through gene expression patterns and co-expression networks, we discovered divergent expression patterns and sub-networks, and found that the members of pair PtrPLATZ14/18 might play different roles in the regulation of macromolecule biosynthesis and modification. Furthermore, we found that PtrPLATZ14 regulates poplar leaf development by affecting cell size control genes PtrGRF/GIF and PtrTCP. In conclusion, our study provides a theoretical foundation for exploring the evolution relationships and functions of the PLATZ gene family within Populus species and provides insights into the function and potential mechanism of PtrPLATZ14 in leaf morphology that were diverse across the Populus genus.

Dynamic nitrogen fixation in an aerobic endophyte of Populus.

Biological nitrogen fixation by microbial diazotrophs can contribute significantly to nitrogen availability in non-nodulating plant species. In this study of molecular mechanisms and gene expression relating to biological nitrogen fixation, the aerobic nitrogen-fixing endophyte Burkholderia vietnamiensis, strain WPB, isolated from Populus trichocarpa served as a model for endophyte-poplar interactions. Nitrogen-fixing activity was observed to be dynamic on nitrogen-free medium with a subset of colonies growing to form robust, raised globular like structures. Secondary ion mass spectrometry (NanoSIMS) confirmed that N-fixation was uneven within the population. A fluorescent transcriptional reporter (GFP) revealed that the nitrogenase subunit nifH is not uniformly expressed across genetically identical colonies of WPB and that only ~11% of the population was actively expressing the nifH gene. Higher nifH gene expression was observed in clustered cells through monitoring individual bacterial cells using single-molecule fluorescence in situ hybridization. Through 15N2 enrichment, we identified key nitrogenous metabolites and proteins synthesized by WPB and employed targeted metabolomics in active and inactive populations. We cocultivated WPB Pnif-GFP with poplar within a RhizoChip, a synthetic soil habitat, which enabled direct imaging of microbial nifH expression within root epidermal cells. We observed that nifH expression is localized to the root elongation zone where the strain forms a unique physical interaction with the root cells. This work employed comprehensive experimentation to identify novel mechanisms regulating both biological nitrogen fixation and beneficial plant-endophyte interactions.

Comprehensive investigation of BZR gene family in four dicots and the function of PtBZR9 and PtBZR12 under drought stress.

Brassinazole-resistant (BZR) transcription factor plays an important role in plant growth and stress resistance through brassinosteroid (BR) signal transduction. However, systematic analysis of the BZR family in dicots remains limited. In this study, we conducted a genome-wide study of four typical dicots: Arabidopsis thaliana, Carica papaya, Vitis vinifera and Populus trichocarpa. Thirty-four BZR gene family members were identified and classified them into three subfamilies. Analysis of promoter and expression patterns revealed crucial role of a pair of homologous BZR genes, PtBZR9 and PtBZR12, in poplar may play a critical role under abiotic stress. PtBZR9 and PtBZR12 were localised in the nucleus and exhibited mutual interactions. Moreover, transient overexpression (OE) of PtBZR9 and PtBZR12 in poplar enhanced tolerance to drought stress. The phenotypic and physiological characteristics of PtBZR9 and PtBZR12 OE in Arabidopsis mirrored those of transient OE in the poplar. Additionally, PtBZR9 and PtBZR12 can bind to the E-box element. Under exogenous BR treatment, transgenic lines displayed a greater decrease in root length than the wild type. Thus, these findings provide a solid foundation for future research on the complex regulatory mechanisms of BZR genes.

Computational dissection of genetic variation modulating the response of multiple photosynthetic phenotypes to the light environment.

BACKGROUND: The expression of biological traits is modulated by genetics as well as the environment, and the level of influence exerted by the latter may vary across characteristics. Photosynthetic traits in plants are complex quantitative traits that are regulated by both endogenous genetic factors and external environmental factors such as light intensity and CO2 concentration. The specific processes impacted occur dynamically and continuously as the growth of plants changes. Although studies have been conducted to explore the genetic regulatory mechanisms of individual photosynthetic traits or to evaluate the effects of certain environmental variables on photosynthetic traits, the systematic impact of environmental variables on the dynamic process of integrated plant growth and development has not been fully elucidated. RESULTS: In this paper, we proposed a research framework to investigate the genetic mechanism of high-dimensional complex photosynthetic traits in response to the light environment at the genome level. We established a set of high-dimensional equations incorporating environmental regulators to integrate functional mapping and dynamic screening of gene‒environment complex systems to elucidate the process and pattern of intrinsic genetic regulatory mechanisms of three types of photosynthetic phenotypes of Populus simonii that varied with light intensity. Furthermore, a network structure was established to elucidate the crosstalk among significant QTLs that regulate photosynthetic phenotypic systems. Additionally, the detection of key QTLs governing the response of multiple phenotypes to the light environment, coupled with the intrinsic differences in genotype expression, provides valuable insights into the regulatory mechanisms that drive the transition of photosynthetic activity and photoprotection in the face of varying light intensity gradients. CONCLUSIONS: This paper offers a comprehensive approach to unraveling the genetic architecture of multidimensional variations in photosynthetic phenotypes, considering the combined impact of integrated environmental factors from multiple perspectives.

Labeled temperate hardwood tree stomatal image datasets from seven taxa of Populus and 17 hardwood species.

Machine learning (ML) algorithms have shown potential in automatically detecting and measuring stomata. However, ML algorithms require substantial data to efficiently train and optimize models, but their potential is restricted by the limited availability and quality of stomatal images. To overcome this obstacle, we have compiled a collection of around 11,000 unique images of temperate broadleaf angiosperm tree leaf stomata from various projects conducted between 2015 and 2022. The dataset includes over 7,000 images of 17 commonly encountered hardwood species, such as oak, maple, ash, elm, and hickory, and over 3,000 images of 55 genotypes from seven Populus taxa. Inner_guard_cell_walls and whole_stomata (stomatal aperture and guard cells) were labeled and had a corresponding YOLO label file that can be converted into other annotation formats. With the use of our dataset, users can (1) employ state-of-the-art machine learning models to identify, count, and quantify leaf stomata; (2) explore the diverse range of stomatal characteristics across different types of hardwood trees; and (3) develop new indices for measuring stomata.

Transcriptome Analysis and Reactive Oxygen Species Detection Suggest Contrasting Molecular Mechanisms in Populus canadensis' Response to Different Formae Speciales of Marssonina brunnea.

Revealing plant-pathogen interactions is important for resistance breeding, but it remains a complex process that presents many challenges. Marssonina leaf spot of poplars (MLSP) is the main disease in poplars; in China, its pathogens consist of two formae speciales, namely, Marssonina brunnea f. sp. Monogermtubi (MO) and M. brunnea f. sp. Multigermtubi (MU). However, the mechanism of the molecular interaction between poplars and the two formae speciales, especially for an incompatible system, remains unclear. In this study, we conducted transcriptome sequencing and reactive oxygen species (ROS) staining based on the interactions between Populus canadensis and the two formae speciales. The results show that the gene expression patterns of P. canadensis induced by MO and MU were significantly different, especially for the genes associated with biotic stress. Furthermore, MO and MU also triggered distinct ROS reactions of P. canadensis, and ROS (mainly H2O2) burst was only observed around the cells penetrated by MU. In conclusion, this study suggested that P. canadensis experienced different resistance reactions in response to the two formae speciales of M. brunnea, providing valuable insights for further understanding the host-pathogen interactions of MLSP.